
*******************************************************
*******************************************************
*******************************************************
*** CREATED BY JONATHAN KING AND JESSICA SCHOENHERR ***
***** REPLICATION DATA FOR "A MATTER OF OPINION?" *****
*******************************************************
*******************************************************
*******************************************************


use "DeathPenaltyDataPart2-20221022.dta"

**********************
*** BASELINE MODEL ***
**********************

reg decisionTherm ib10.treatmentGroup
est store dpBaseline

est restore dpBaseline
margins, at(treatmentGroup=(6(1)10))

est restore dpBaseline
margins, at(treatmentGroup=(10)) ///
	at(treatmentGroup=(6)) post
test _b[1._at] = _b[2._at]
* stat sig (0.0303)

est restore dpBaseline
margins, at(treatmentGroup=(10)) ///
	at(treatmentGroup=(7)) post
test _b[1._at] = _b[2._at]
* not sig

est restore dpBaseline
margins, at(treatmentGroup=(10)) ///
	at(treatmentGroup=(8)) post
test _b[1._at] = _b[2._at]
* not stat sig

est restore dpBaseline
margins, at(treatmentGroup=(10)) ///
	at(treatmentGroup=(9)) post
test _b[1._at] = _b[2._at]
* not stat sig

est restore dpBaseline
margins, at(treatmentGroup=(6)) ///
	at(treatmentGroup=(7)) post
test _b[1._at] = _b[2._at]
* not stat sig

est restore dpBaseline
margins, at(treatmentGroup=(6)) ///
	at(treatmentGroup=(8)) post
test _b[1._at] = _b[2._at]
* not stat sig

est restore dpBaseline
margins, at(treatmentGroup=(6)) ///
	at(treatmentGroup=(9)) post
test _b[1._at] = _b[2._at]
* stat sig = 0.0072

est restore dpBaseline
margins, at(treatmentGroup=(7)) ///
	at(treatmentGroup=(8)) post
test _b[1._at] = _b[2._at]
* not stat sig

est restore dpBaseline
margins, at(treatmentGroup=(7)) ///
	at(treatmentGroup=(9)) post
test _b[1._at] = _b[2._at]
* stat sig = 0.0522

est restore dpBaseline
margins, at(treatmentGroup=(8)) ///
	at(treatmentGroup=(9)) post
test _b[1._at] = _b[2._at]
* stat sig = 0.0433

********************
*** PARTISANSHIP ***
********************

* dem 1
* ind 2
* rep 3

reg decisionTherm ib10.treatmentGroup##ib1.partisanship
est store dpPartisanship

est restore dpPartisanship
margins, at(partisanship=(1(1)3) treatmentGroup=(6(1)10))

* democrats across treatment
*compared to group 6
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0791

est restore dpPartisanship
margins, at(treatmentGroup=(7) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 7
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(7) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(1)) ///
	at(treatmentGroup=(7) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0685

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(1)) ///
	at(treatmentGroup=(7) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0484

*compared to group 8
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(8) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(1)) ///
	at(treatmentGroup=(8) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 9
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(9) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

* republicans across treatment
* compared to group 6
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(3)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0200

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(3)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0073

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(3)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(7) partisanship=(3)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 7
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(3)) ///
	at(treatmentGroup=(7) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(3)) ///
	at(treatmentGroup=(7) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(3)) ///
	at(treatmentGroup=(7) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 8
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(3)) ///
	at(treatmentGroup=(8) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0121

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(3)) ///
	at(treatmentGroup=(8) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0039

*compared to group 9
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(3)) ///
	at(treatmentGroup=(9) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

* dems vs. reps within a treatment
est restore dpPartisanship
margins, at(treatmentGroup=(6) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(7) partisanship=(1)) ///
	at(treatmentGroup=(7) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0033

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(1)) ///
	at(treatmentGroup=(8) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(1)) ///
	at(treatmentGroup=(9) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0046

est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(10) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0013

**************
*** GENDER ***
**************

reg decisionTherm ib10.treatmentGroup##i.female
est store dpGender

est restore dpGender
margins, at(treatmentGroup=(6(1)10) female=(0(1)1))

* men across treatment
* compared to group 1
est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(6) female=(0)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0031

est restore dpGender
margins, at(treatmentGroup=(9) female=(0)) ///
	at(treatmentGroup=(6) female=(0)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0117

est restore dpGender
margins, at(treatmentGroup=(8) female=(0)) ///
	at(treatmentGroup=(6) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(7) female=(0)) ///
	at(treatmentGroup=(6) female=(0)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 7
est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(7) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(0)) ///
	at(treatmentGroup=(7) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(8) female=(0)) ///
	at(treatmentGroup=(7) female=(0)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 8
est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(8) female=(0)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0522

est restore dpGender
margins, at(treatmentGroup=(9) female=(0)) ///
	at(treatmentGroup=(8) female=(0)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 9
est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(9) female=(0)) post
test _b[1._at] = _b[2._at]
* no

* women across treatment
*compared to group 6
est restore dpGender
margins, at(treatmentGroup=(10) female=(1)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(1)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(8) female=(1)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(7) female=(1)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 7
est restore dpGender
margins, at(treatmentGroup=(10) female=(1)) ///
	at(treatmentGroup=(7) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(1)) ///
	at(treatmentGroup=(7) female=(1)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0756

est restore dpGender
margins, at(treatmentGroup=(8) female=(1)) ///
	at(treatmentGroup=(7) female=(1)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 8
est restore dpGender
margins, at(treatmentGroup=(10) female=(1)) ///
	at(treatmentGroup=(8) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(1)) ///
	at(treatmentGroup=(8) female=(1)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 9
est restore dpGender
margins, at(treatmentGroup=(10) female=(1)) ///
	at(treatmentGroup=(9) female=(1)) post
test _b[1._at] = _b[2._at]
* no

* men vs women, within a treatment
est restore dpGender
margins, at(treatmentGroup=(6) female=(0)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no 

est restore dpGender
margins, at(treatmentGroup=(7) female=(0)) ///
	at(treatmentGroup=(7) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0350

est restore dpGender
margins, at(treatmentGroup=(8) female=(0)) ///
	at(treatmentGroup=(8) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(0)) ///
	at(treatmentGroup=(9) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(10) female=(1)) post
test _b[1._at] = _b[2._at]
* stat sig 0.0015

******************
*** FULL MODEL ***
******************

reg decisionTherm ib10.treatmentGroup##ib1.partisanship##i.female
est store dpFull

est restore dpFull
margins, at(treatmentGroup=(6(1)10) partisanship=(1(1)3) female=(0(1)1))

* democrat men across treatment
*compared to group 6
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* not

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0732, 63.8 vs. 55.7

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

* democrat women across treatment
*compared to group 6
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no 

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

* Republican men across treatment
*compared to group 6
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0162

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes; p = 0.0087

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0647

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes; p = 0.0409

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0089

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*yes; p = 0.0042

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

* Republican women across treatment
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*no

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

* democrat men vs women, within a treatment
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(6) partisanship=(1))

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(7) partisanship=(1))

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(8) partisanship=(1))

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(9) partisanship=(1))

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0165

est restore dpFull
margins, dydx(female) at(treatmentGroup=(10) partisanship=(1))

* republican men vs women, within a treatment
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(6) partisanship=(3))

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(7) partisanship=(3))

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(8) partisanship=(3))

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(9) partisanship=(3))

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(10) partisanship=(3))

* men, republican vs. democrat, within a treatment
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(6) female=(0))

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(7) female=(0))

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no, p = 0.0754

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(8) female=(0))

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0010

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(9) female=(0))

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0001

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(10) female=(0))

* women, republican vs. democrat, within a treatment
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(6) female=(1))

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(7) female=(1))

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(8) female=(1))

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, 0.0491

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(9) female=(1))

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(10) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, 0.0051

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(10) female=(1))

***************************
*** WITH ADDED CONTROLS ***
***************************

* destring income, education, and age
replace income = "." if income == "NA"
destring(income), replace

replace education = "." if education == "NA"
destring(education), replace

reg decisionTherm ib10.treatmentGroup##ib1.partisanship##i.female income education age
est store dpFullControl

est restore dpFullControl
margins, at(treatmentGroup=(6(1)10) partisanship=(3) female=(0(1)1))

est restore dpFullControl
margins, dydx(female) at(treatmentGroup=(10) partisanship=(3))

est restore dpFullControl
margins, dydx(partisanship) at(treatmentGroup=(10) female=(1))

* democrat men across treatment
*compared to group 6
est restore dpFullControl
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFullControl
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFullControl
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFullControl
margins, at(treatmentGroup=(7) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFullControl
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* not

est restore dpFullControl
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFullControl
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 8
est restore dpFullControl
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFullControl
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 9
est restore dpFullControl
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

* democrat women across treatment
*compared to group 6
est restore dpFullControl
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no 

est restore dpFullControl
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFullControl
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFullControl
margins, at(treatmentGroup=(7) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFullControl
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFullControl
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFullControl
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 8
est restore dpFullControl
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFullControl
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 9
est restore dpFullControl
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

* independents men across treatment
*compared to group 6
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
* not

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0732, 63.8 vs. 55.7

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]
* no

* independent women across treatment
*compared to group 6
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no 

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no

* independents across gender
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(2) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

* independent men vs. dem men
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

* rep men vs. ind men
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(2) female=(0)) post
test _b[1._at] = _b[2._at]


* dem women vs. independent women
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(10) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

* rep women vs. independent women
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(10) partisanship=(2) female=(1)) post
test _b[1._at] = _b[2._at]


